Flow equations in the light-front QCD: mass gap and confinement 
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The light-front QCD is studied using the method of flow equations. Solving 
the light-front gluon gap equation, the effective gluon mass is generated dynam- 
ically. The effective interaction between static quark and antiquark, generated 
through elimination of the quark-gluon minimal coupling by flow equations, has 
the Coulomb, 1/q 2 , and confining, 1/q 4 , singular behavior. Elimination of the 
quark-gluon coupling at small gluon momenta is governed by the cutoff depen- 
dent, dynamical gluon mass, which makes this elimination possible and provides 
such an enhancement at q ~ 0. The cutoff, which regulates small light-front x 
divergences, sets up a scale for the dynamical gluon mass and the string tension in 
the gq-potential. The mechanism of confimenemet in the light-front frame is sug- 
gested, based on the singular nature of the light-front gauge along the light-front 



1 Introduction 

The perturbative aspects of non-abelian gauge theories were underestood many 
years ago, and the perturbative calculations provided convincing proof that 
QCD is the theory of strong interactions. However nonperturbative QCD phe- 
nomena have been difficult to analyze mainly because calculational techniques 
are still lacking, even though the qualitative features have been more or less 
understood. 

In particular, it is widely believed that pure Yang- Mills theory, with no dy- 
namical quarks, posseses the features like asymptotic freedom, mass generation 
through the transmutation of dimensions and confinement: linear potential be- 
tween static (probe) quarks. In the past few years there were several studies 
addressing the issue of confinement and generation of mass gap in the frame- 
work of the Schrodinger picture , 2 . To mention a few early works, see refs. 
3 . We have tried to understand these nonperturbative mechanisms also in a 
Hamiltonian framework, solving flow equations for canonical QCD Hamilto- 
nian in the light-front (LF) gauge selfconsistently for the few lowest sectors. 
Dynamics of quarks has been excluded to disentangle the complexity of chiral 
symmetry breaking. Early studies of confinement in the LF frame were con- 
ducted using similarity renormalization 4 and transverse lattice 5 , based on the 
fact that QCD in 3 + 1 already has a confining interaction in the form of the 
instantaneous four-fcrmion interaction, which is the confining interaction in 
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1 + 1. However, in our study the instantaneous interaction is canceled by the 
corresponding term in the dynamical interaction, generated by flow equations. 
The result is a three dimensional linear rising potential in the infrared region. 

Our basic strategy has been to generate dynamically an effective gluon 
mass (energy) through interaction with the LF 'vacuum' (exact zero modes on 
the light-front, x = 0), and construct an effective interaction between static 
quarks as an exchange with this dynamical gluon, which forms a flux between 
quarks. One could construct a LF field theory with massive photons and glu- 
ons, as was done in 6 , but this theory is in a Coulomb phase rather than confin- 
ing because of the lack of a dynamical mass generation mechanism. Dynamical 
mass generation in the light-front frame dates back to the works of Cornwall 
7 ', where a kind of pinch technique was used to define a gluon propagator and 
hence a gluon effective energy, extracted from a Green's function for physical 
observables (for example two-body scattering amplitude). It was suggested in 
7 , that in physical terms, a gluon 'mass' may lead to a vortex condensation a 
with < TrG 2 ^ >^ 0, and conversely < TrG 2 ^ >^ generates a gluon mass. 
This provides a link with the existing vortex picture of confinement. 

A subtle point in the LF field theory is that the LF vacuum is just empty 
space. Therefore it seems a problem how confinement can occur in the LF 
frame, and what quantity sets up a scale for a dynamical mass gap and the 
string tension. We adopt a model, suggested by Susskind and Burkardt in the 
context of chiral symmetry breaking in the LF frame 8 . In the parton model 
one pictures a fast moving hadron as being some collection of constituents with 
relatively large momentum, such that when one boosts the system, doubles its 
momentum, all these partons double their momenta and so forth. Therefore 
one can formulate a field theory on the axis of the light-front momentum x. 
Partons that form a hadron are at positive, finite x and according to Feynman 
and Bjorken fill the x-axis in a way which gets denser and denser as one goes 
to smaller x; and the vacuum is at x — 0. The fundamental property of LF 
Hamiltonians, that under a rescaling of the LF momentum, x — > Ax, the LF 
Hamiltonian scales like H — ► H/ A, can be interpreted as a dilatation symmetry 
along the x-axis, if one thinks of the x-axis as a spacial axis. This symmetry 
holds on a classical level and it is brocken on a quantum level by anomali. As 
one approaches small x, interaction between partons gets stronger, contributing 
divergent matrix elements. A natural cutoff is provided by Sx = ex, a minimal 
spacing between constituents, which plays the role of UV-regulator b . As long 
as Sx is finite, i.e. as long as the density of partons on the x-axis is not infinite, 

a In the massive gauge theory one needs to introduce a scalar field, which condenses, in order 
to preserve gauge invariance. Topologically this field represents vortices 7 . 
''Small x correspond to the large light-front energies. 
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one obtains finite matrix elements. Cutoff 6x breaks the dilatation symmetry 
along the x-axis and gauge symmetry, and sets up a scale for quantities of 
dimension of energy (mass) and higher powers in energy. In particular, one 
generates dynamically an effective gluon mass, which in turn will define string 
tension between quark and antiquark. Formation of the qq bound state through 
breaking an internal symmetry can be viewed analogously to the creation of 
Cooper pairs in superconducter. 

In terms of effective theory, the Hamiltonian below the LF cutoff, H< e , 
which desribes high energy (UV) strongly correlated states, can be substi- 
tuted by the v.e.v., since strongly coupled configurations are frozen. Hamilto- 
nian above the cutoff H >e is treated by flow equations (renormalization group 
transformation). One way of looking at the physics behind this v.e.v. and 
mass generation in gluodynamics is that the composite field <j> which creates 
+ glueballs has a finite v.e.v. 7 , 2 . 

The dilatation symmetry, discussed above, reflects some underlying scale 
invariance of the LF field theory formulated on x-axis. Introducing the cutoff, 
breaks this symmetry. Because physics should remain independent on the 
cutoff, one must be looking for a fixed point of the renormalization group. 
Therefore the right tool for studing such a system is the renormalization group, 
which is provided by the method of flow equations 9 . 

In the method of flow equations an effective interaction between quarks 
arises through elimination of quark-gluon coupling. Procedure converges in the 
UV for large gluon momenta transfer q 2 , but in the IR one is not able to elim- 
inate vanishing gluon momenta, because the similarity factor which governes 
the elimination does not decay for vanishing arguments. In physical terms, one 
can not integrate soft gluons in the same fashion. As soon as gluon aquircs a 
dynamical mass, which vanishes at large gluon momenta and not equal zero 
only at small momenta, the similarity factor decays even for small momenta 
transfer, because the argument contains an effective gluon energy instead of 
only gluon momctum. Therefore we can eliminate by flow equations even soft 
gluon modes, that are responsible for the long-range part of an effective quark 
potential. In the IR this gives 1/q 4 behavior for an effective qq interaction. 

In the next section we solve flow equations for an effective gluon mass and 
quark-antiquark effective interaction, based on the QCD Hamiltonian in the 
light-front gauge. 

2 Gluon mass gap and an effective quark interaction 

Let Q being a projection operator on a one-gluon state, and P on a qq state, 
Q\ip) = \g) and P\ip) = \qq). Flow equations for matrix elements of the QCD 
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Hamiltonian between these states read 



dh pp ,{l) _ \ ^ / dh pq (l) 1 1 dh qpl (l) 



dl ^\ dl E p (l)-E q (l) hqp ' {l) + hpqil) E p ,(l)-E q (l) dl 

where p (p') runs through all states in the P-subspace, and q in the Q-subspace. 
Flow equations for the quark-gluon coupling h pq (coupling between P and Q 
sectors) are 



dh pq (l) 
dl 



(E p (l)-E q (l)) Vpq (l) 
h pq {l) d 



z pq (l) = l(E p (l)-E q (l)f , (2) 

where / is the similarity function with properties /(0) = 1 and f(z — ► oo) = 0. 
From Eq. (2), elimination of quark-gluon coupling 

h pg (l) = h pq (0)f(z pg (l)), (3) 

is governed by the similarity function f(l(E p (l) — E q (l)) 2 ), which vanishes for 
the matrix elements with energy differences exceeding the cutoff A, \E p (l) — 
E q (l)\ >• \j\fl = A, h pq (0) is the initial value, and E q (l) (E p (l)) is a solution of 
the flow (renormalization group) equation, Eq. (1). As long as the argument 
of the similarity function is not equal zero, one can eliminate quark-gluon 
coupling and solve flow equations for an effective gluon energy and effective 
qq-interaction, Eq. (1). However, if one neglects the cutoff dependence of 
energies, E p (l) <~ E p (0), the argument of the similarity function may become 
zero in the degenerate case, E p (0) = E q (0), and then the energy denominator 
blows up in effective interactions in Eq. (1). 

In physical terms, the leading behavior of the argument in the quark-gluon 
similarity function is given by an effective gluon energy transfer between quark 
and antiquark, (E p (l) — E q (l)) ~ E e J^ on (l), which reduces at the initial value 

of flow parameter, I = 0, to a gluon momentum transfer, E e gl ^ on (Q) ~ q 2 . 
Neglecting the cutoff dependence of the effective gluon energy, one eliminates 
quark-gluon coupling only at large gluon momenta, and fails at small gluon 
momenta, because the argument of the similarity function tends to zero. For 
small gluon momenta the dependence of gluon effective energy on the cutoff 



4 



becomes important, given by a dynamically generated gluon mass, that makes 
elimination of quark-gluon coupling possible even for small gluon momenta. 
We therefore solve flow equations, Eq. (1), consistently for an effective gluon 
mass and an effective gg-interaction below. We show, that elimination of quark- 
gluon coupling for high and low gluon momenta contribute correspondingly to 
the UV and IR parts of quark-antiquark potential, that behave in momentum 
space like 1/q 2 and 1/q 4 , respectively, where q is a gluon momentum transfer. 



2.1 Gluon gap equation 

In the light-front frame flow equation for an effective gluon mass (the first 
equation in Eq. (1)), with connection to the light-front energy q~ = {q\ + 
fi 2 (\))/q + , reads 

<^ 2 (A) orT1 Ar f 1 dx f°° d 2 k ± 2/xN 1 df 2 {Ql{\);\) 



f L dx f°° d 2 k ± 



d\ 1 J Jo x(l-x)J 16tt 3 ^ v Q\{X) d\ 

kl+m 2 



x(l — x) 



Ww 3 " 9 " 'Qi(X) dX 
CI - xf 



x + l + — i— )) , (4) 



with 



where gluon couples to the quark- anti-quark pairs, with the quark-gluon cou- 
pling g q (X), and pairs of gluons, with the three-gluon coupling g g (X); for finite 
A these couplings are different from each other functions of momenta; con- 
nection between flow parameter and the cutoff, I = 1/A 2 , is used; similarity 
function / plays the role of UV regulator in loop integrals; current quark 
mass is m; for SU(N C ) T f 5 ab = Tr(T a T b ) = \8 ab and the adjoint Casimir is 
C a Sab = f acd f bcd = N c S abl N c is the number of colors (i.e., N c = 3). 

Coupled system of equations for the cutoff dependent, effective gluon and 
quark masses was first derived in 11 . Here, neglecting the cutoff dependence of 
quark mass this system is reduced to Eq. (4). 

Generally, it is difficult to solve Eq. (4) . One of the reasons is that this 
equation contains (unknown) running coulings; therefore the cutoff dependence 
of couplings is neglected below. Second, initial condition for Eq. (4) is not 
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known. The following renormalization condition to define the running gluon 
mass /i(A) through the 'physical' mass is assumed 11 : the effective Hamiltonian 
at the scale A has bosonic eigenstates with eigenvalues of the form q~ = (q 2 ± + 
AW", i-e. 

Q± + M 2 / / \ gj + M 2 (A) ; , . 
—^{q,q) = ^ (q,Q) 

- f dX' J2 (vMWhvW - h qlp (\') Vpq (\')) , (6) 
J p 

where jl denotes the 'physical' gluon mass; the generator rj pq , given by Eq. (2), 
eliminates the quark-gluon (three-gluon) coupling h pq ; \q) denotes a single 
effective gluon state with momentum <7 + and q_\_, (q'\q) = I6ir 3 q + S^ (q' — q). 
Neglecting the cutoff dependence of the gluon mass on the r.h.s. of Eq. (4), 
the renormalization point is choisen at q 2 = jl 2 (for in- and out-going and 
intermediate state gluons). Third, severe divergences arise at small light-front 
momenta x and must be regularized. For this purpose, as suggested by Zhang 
and Harindranath in 10 , the minimum cutoff u for transverse momentum k± 
is introduced. Therefore, flow equation, Eq. (4), is integrated in the finite 
limits, [it; A]. Explicitly, from the similarity function in gluon loop (the second 
term in Eq. (4)) one has u < Q\ = {k\ + jl 2 )/x(l - x) - fi 2 < A, that 
restricts the transverse momentum to k± m i n = (u 2 + fi 2 )x{l — x) — ft 2 < 
k± < k± max = (A 2 + fi 2 )x(l — x) — jl 2 , and the light-front momentum to 
p. 2 /(u 2 + p?) < x < 1 — fi 2 /(u 2 + jl 2 ). In the same way one finds restrictions 
on mometa in the quark-gluon loop (the first term in Eq. (4)). Integrating 
Eq. (4) with all constraints discussed above and assuming the condition of 
Eq. (6), gives 12 

ji 2 (A) = fl 2 + 5n 2 PT (A) + 5n 2 dyn (A, A, u) . (7) 
where the perturbative term 

^ T W = £,X 2 (c a Onf 2 -^) + T f N f ^j . (8) 

reproduces the result of the LF perturbation theory, with renormalization 
point q 2 — 0, 10 . When renormalization is performed in the second oder, the 
perturbative mass correction is absorbed by the mass counterterm, mr CT = 
— <5^p T (A = A — > oo). After the perturbative renormalization is completed, 
the dynamical mass 



A 2 

l4yn( X ) = M 2 + fydyn( X , M> u ) = + u ) ln ~ 
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») = -iJf 2 + T,N I{ \ + fi). (9) 

is left. Note, that the gluon mass [idyn is generated dynamically by flow equa- 
tions. In the limit /} — > 0, one has er = lim^ it) = u 2 g 2 C a /2TT 2 , and, 
as shown below, a plays the role of the string tension between quark and 
antiquark. Scale u is introduced by regulating the divergences at small LF 
momenta, x <~ 0, (small LF a; rcgularization). The value u sets up a scale for 
the dynamical gluon mass and the string tension. 



2.2 Effective qq-interaction 

An effective interaction, generated between quark and antiquark by flow equa- 
tions in the LF frame (the second equation in Eq. (1)), reads 

V qS = -Const^Y) lim B„ v , (10) 

ji— >o 

where instead of the strong coupling a s some constant, Const, is introduced; 
( r y^j v ) is a current-current term in the exchange channel. At the end we set the 
gluon mass renormalization point ('physical' gluon mass) to zero, p, —> 0; thus 
eliminating the sensitivity to the renormalization point. The kernel includes 
dynamical interaction and instantaneous exchange with nonzero dynamical 
gluon mass, given by 

SQ 2 

B^ v = g^ v (h + h) + VnVv—p^ (h - h) , (11) 

where g^ v is the LF metric tensor, and the LF vector 77^ is defined as r\- k = k + . 
The cutoff dependence of four-momentum transfers along quark and antiquark 
lines is accumulated in the factor 

with mometum transfers defined in the light-front frame as 

QiW = Ql + mUa) , Qi = ^-^ + ^-*') 2 

y XX 
Q\(\) = Ql+ fJ-dynW ' Ql = QlU^(l-x); s'->(l-x') 7 (13) 

and the dynamical gluon mass ^dyn is given by Eq. (9); also the following 
momenta are used Q 2 = (Ql + Ql)/2 and 5Q 2 = (Ql — Q%)/2. Calculating 
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Eq. (10) with an explicit form of similarity function, gives in the leading order 
5Q 2 <C Q 2 the following gg-interaction 12 



where to this order Q 2 reduces in the instant frame to the square of gluon mo- 
mentum transfer, Q 2 ~ q 2 , with the gluon momentum q = (q z ,q±). In order 
to reproduce the standard Coulomb and linear confining potentials, the correct 
prcfactors are restored, using the freedom to fit Const and a terms. Confining 
term in Eq. (14), with singular behavior like l/<f 4 , arise from the elimination 
of the quark-gluon coupling at small gluon momenta, that is governed by the 
cutoff dependent, dynamical gluon mass. 
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